home *** CD-ROM | disk | FTP | other *** search
-
-
-
- SSSSSSSSIIIIDDDDIIII((((3333FFFF)))) SSSSSSSSIIIIDDDDIIII((((3333FFFF))))
-
-
-
- NNNNAAAAMMMMEEEE
- SSIDI - SSIDI computes the determinant, inertia and inverse of a real
- symmetric matrix using the factors from SSIFA.
-
-
- SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
- SUBROUTINE SSIDI(A,LDA,N,KPVT,DET,INERT,WORK,JOB)
-
- DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
- On Entry
-
- AAAA REAL(LDA,N)
- the output from SSIFA.
-
- LLLLDDDDAAAA INTEGER
- the leading dimension of the array A.
-
- NNNN INTEGER
- the order of the matrix A.
-
- KKKKPPPPVVVVTTTT INTEGER(N)
- the pivot vector from SSIFA.
-
- WWWWOOOORRRRKKKK REAL(N)
- work vector. Contents destroyed.
-
- JJJJOOOOBBBB INTEGER
- JOB has the decimal expansion ABC where
- If C .NE. 0, the inverse is computed,
- If B .NE. 0, the determinant is computed,
- If A .NE. 0, the inertia is computed.
- For example, JOB = 111 gives all three. On Return Variables not
- requested by JOB are not used.
-
- AAAA contains the upper triangle of the inverse of
- the original matrix. The strict lower triangle
- is never referenced.
-
- DDDDEEEETTTT REAL(2)
- determinant of original matrix.
- Determinant = DET(1) * 10.0**DET(2)
- with 1.0 .LE. ABS(DET(1)) .LT. 10.0
- or DET(1) = 0.0.
-
- IIIINNNNEEEERRRRTTTT INTEGER(3)
- the inertia of the original matrix.
- INERT(1) = number of positive eigenvalues.
- INERT(2) = number of negative eigenvalues.
- INERT(3) = number of zero eigenvalues. Error Condition
-
- AAAA division by zero may occur if the inverse is requested and SSICO has
- set RCOND .EQ. 0.0 or SSIFA has set INFO .NE. 0 . LINPACK. This
-
-
-
- PPPPaaaaggggeeee 1111
-
-
-
-
-
-
- SSSSSSSSIIIIDDDDIIII((((3333FFFF)))) SSSSSSSSIIIIDDDDIIII((((3333FFFF))))
-
-
-
- version dated 08/14/78 . James Bunch, Univ. Calif. San Diego, Argonne
- Nat. Lab Subroutines and Functions BLAS SAXPY,SCOPY,SDOT,SSWAP Fortran
- ABS,IABS,MOD
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- PPPPaaaaggggeeee 2222
-
-
-
-